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Abstract. We study a one-dimensional system composed of three charged bosons confined in an external 
harmonic potential. More precisely, we investigate the ground-state correlation properties of the system, 
paying particular attention to the strong- interaction limit. We explain for the first time the nature of 
the degeneracies appearing in this limit in the spectrum of the reduced density matrix. An explicit rep- 
resentation of the asymptotic natural orbitals and their occupancies is given in terms of some integral 
equations. 
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1 Introduction 

In recent years there has been a growing interest in sys- 
tems of interacting particles trapped in potential wells be- 
cause of their possible use in quantum information tech- 
nology [T] . Especially, systems composed of particles held 
together in harmonic potentials have drawn considerable 
theoretical attention |2l3l4l5l6l7l8l9ll0lllll2ll3ll4ll5ll6ll71 
|18I19I20I21I22I23I24I25I26I27I28J . Besides fermionic sys- 
tems, which serve well as models of quantum dots, many 
attempts have been made to explore the properties of the 
bosonic ones with a contact potential |15I16I17I18I19I20|2T] . 
Such systems in the one-dimensional (ID) limit have at- 
tracted much attention and their properties appear to 
be well understood |17ll8ll9l20l2T] . Also, there has been 
considerable interest in the properties of artificial atoms 
composed of bosons interacting via a Coulomb potential 
|22I23I24I25I28I26I27J . which in turn serve well as models 
of electromagnetically trapped ions [55]. However, there 
has been so far relatively little theoretical research on such 
systems in the strictly ID limit. 

In the present paper we consider an ideal ID sys- 
tem composed of three identical bosons described by the 
Hamiltonian 

H = 



where Xij = Xi — x 
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spatial variable 



J are the interparticle distances and 
2:3) is the centre-of-mass. In ([T|), the 



X is given in oscillatory units ^/K/rrwj, 
and g is the ratio of the Coulomb and the confinement 

energies: g = 4y^- 

Experimentally, a ID configuration of trapped ions 
can be realized in a 3D harmonic trap with a transverse 
trapping frequency a;_L much larger than the axial one 
w, e = ijj±/ijj ^ 1. Although the ID Hamiltonian ((T)) is 
strictly valid only in the limit e ^- 00 , it works well even at 
moderate values of e if the confinement in the x direction 
is very weak {g — >■ 00). For more details on this point we 
refer the reader to [H]. 

The main goal of this paper is to make a detailed inves- 
tigation of the ground-state correlation properties of the 
system ([T|) in the strong- interaction limit {g — > 00). To this 
end, we use an approximation based on the second-order 
Taylor series expansion of ^ around the classical equilib- 
rium distances of the particles |28| . Within the framework 
of this approximation, an asymptotically exact expression 
for the ground-state bosonic wavefunction at the 5 — > 00 
limit can be obtained. We derive an explicit expression 
for the asymptotic reduced density matrix (RDM) and in- 
vestigate for the first time the nature of the degeneracies 
appearing in its spectrum. We provide an explicit repre- 
sentation of the asymptotic natural orbitals and their oc- 
cupancies, given by integral equations independent of g. 
In particular, we show that only the three natural orbitals 
contribute significantly to the asymptotic bosonic ground- 
state. Moreover, to gain insight into the general features 
of the Schmidt expansion of the RDM, we determine nu- 
merically the values of the three lowest occupancies over 
a wide range of values of g. 

This paper is organized as follows. Section [5] derives 
closed-form analytical approximate solutions for the bosonic 
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ground-state of the system ([TJ. Section |3] tests their vahd- 
ity and provides, in particular, detailed results for the de- 
pendencies of the degree of correlation on g. Section U] de- 
rives g-indcpendcnt integral equations defining the asymp- 
totic natural orbitals and their occupancies. Here, two dif- 
ferent forms of the Schmidt expansion of the asymptotic 
RDM will be discussed. Finally, some concluding remarks 
are placed in Section [5] 



2 Harmonic approximation 

The potential ^ attains its minimum at six points, namely 
(x™,x™,a;™), X™ = — Xc, x™ = 0, x™ ~ Xc, where Xc ~ 

(105)3/2, and {i,j,k} are the permutations of {1,2,3}. 
Consequently, the distances between the three particles in 
the ground-state of the system in the classical limit are 



so that the corresponding Schrodinger equation takes the 
form 

^^ 2dxy '""^ 2^ 3 > + 
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+ — (is - Xl)^}lpap = Eap-)pap- (10) 
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Eq. (fTO|) would suggest an ansatz of the form 



e *• 3 -' e 



i'ap{S:i,X2,X3) = 



a;^ = (10g)3/2, cc™ = x^ - a;™ 



32 ~ -^3 -^2 ~ 

^ — 2xc = (103)3, and 



X™ = (x™ -I- X™ -I- 2:3" )/3 = 0, where we have referred to 
the minimum {i = 1, j = 2,k = 3}. The Taylor expansion 
of © around (a;^J,a;^2: a^3i.-^") is 



i>]=l 



dv ^ 

-g^ Uij=xg (a^ij - X,j) + 



2^2,2 la:>j=a:J5 (a;jj X^^ ) +...+ 



Noticing that 



dV 



7 . a^.. Ui3=a;j;j (a^i a;^- a;^) — 0, 
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(■jp^ L -2;m= (— l)*^"'-^^-T-) and retaining the terms up 
to second order, one gets the approximation 

29 29 

+^(a;2 -Xi~ Xcf + —{X3 -X2- Xcf + 
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Thus the Hamiltonian ([1} is approximated by 
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It is convenient to introduce new coordinates in ([5]), 

Xi l-^ Xi - Xc,X2 1-^ X2,X3 K^ £'3 +Xc, (9) 



After substituting (|lip into pop . and performing some 
straightforward algebra, we find parameter values a, /?, 7 
at which the function pT|) satisfies Eq. PH)) . that is 



^^4vf'^=|j^^^^-^)' 



and 



^ap = ^(10 + 10\/3 + 2%/l45+ 15* 103gi). (12) 

Changing the variables back in ()lip in accordance with 
(|ni), we can construct then the approximate spatial sym- 
metric wavefunction: 

-0b''(xi,X2,X3) = C ^ '(/'ap(a;i -|-Xc,Xj,Xfc - Xc),(13) 

where {i,j, fc} are again the integers {1,2,3} permuted 
into a different order. The normalization constant in ([T5| 
can be calculated analytically. However, since it is a quite 
lengthy formula, we report here only its value as (7 —^ 00 



a 



(f)^ 



-\/238 7r5 



(14) 



Before going further it should be stressed that the ap- 
proximation ([5]) coincides with that obtained from the 
second-order Taylor series expansion of ^ around the 
classical equilibrium positions of the particles (— Xc, 0, Xc). 
The approximation strategy considered here yields thus 
the results consistent with those of the standard normal- 
mode theory [2S] . 



3 Numerical tests 

The system of Bose particles described by the Hamiltonian 
^ gets fermionized for any 5 7^ J21I27| and its ground- 
state wavefunction V's can be related to the lowest energy 
antisymmetric wavefunction ipp by 



V'b(xi,X2,X3) = |V'_f(xi,X2,X3)|. 



(15) 
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Therefore, in the hmit g -^ 0, V'b tends to the modulus of 
the Slater determinant 



1 



V"^ {xi,X2,X3) = -=\det„^Q.j^^{ip„{xj))\, (16) 



where {(pn} are the single-particle orbitals of the ideal sys- 
tem {g = 0). In contrast, in the non- interacting case, ipB is 
nothing else but the product function ^^ = Y[i=i foi^i)- 
The above in turn implies that the quantities associated 
with ipB exhibit generally a discontinuity at the point 
g = 0. In particular, when it comes to the energy, it tends 
to 4.5 as g — >■ 0, while at g = 0, it has the value 1.5. In or- 
der to test the applicability of the approximations P^ and 
(fT3|) . we determined numerically the ground-state bosonic 
wavefunction ipB and its corresponding energy, for a wide 
range of values of g. The results of Eq. (fT^ arc compared 
with our accurate numerical results in Fig. [U from which 
it can be seen that the approximate energy tends from be- 
low to the exact one with increasing g. Surprisingly, Eq. 
(fT2|) yields very good estimates of the true values already 
when g exceeds the value g ~ 5. However, the results of 
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Fig. 1. Comparison of the approximations ([T^ (bro- 
ken curve) with the numerically exact results (continuous 
curve) . 



Fig. [T] are not a good indicator of the accuracy of the 
approximate bosonic wavefunctions (|13p. To gain insight 
into their range of applicability, we analyse their ability to 
reproduce the numerically exact bechaviour of the degree 
of correlation 1301 



K^[trp']-'^[Jp{x,yf 



dxdy] 



where p is the RDM expressed in coordinates 



p{x,y)= / 'il;{x,X2,X3)'i{;{y,X2,X3)dx2dx3. 



(17) 



(18) 



The degree of correlation K counts approximately the 
number of orbitals actively involved in the Schmidt de- 
composition of the RDM and is one of the transparent 
measures of correlation. It is worth stressing at this point 



that the linear correlation entropy L, which is also a pop- 
ular measure of correlation |6l7l8llll32] . is related to K 
via L = 1 — 1/K. The results for the degree of correla- 
tion calculated from the numerically exact bosonic wave- 
functions and the approximate ones -0^^ p3p are plot- 
ted in Fig. [5] up to g = 200. As one can see, accept- 
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Fig. 2. Degree of correlations calculated from the nu- 
merically exact wavefunctions (continuous curve) and the 
approximate ones (|13p (broken curve), as a function of 
Ing. The horizonal line marks the result determined as 
g ^- with the use of dH 



able results are reached just at a value of about g — 200 
(In 200 sa 5.29). To complete our discussion, we compare 
the densities n{x) = -ip^x, x) evaluated from the approxi- 
mate bosonic wavefunctions (J13p with the exact ones de- 
termined numerically. Our results are displayed in Fig. 




Fig. 3. Comparison of one-particle density n{x) = 
3p{x, x) calculated from the approximate bosonic wave- 
function ([T^ (broken curve) with the exact one deter- 
mined numerically (continuous curve), for g — 20 (In 20 ~ 
3), 5 = 50 (In 50 w 3.9), and g = 200 (In 200 = 5.29). 



121 for three different values of the interaction strength: 
g = 20, 50, and g = 200. As one could have expected, the 
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agreement between the approximate and exaet densities equilibrium positions of the particles. It is worth empha- 
is observed, at least within the graphical accuracy, only in sizing that the peaks at x = ±Xc have the same profile, 
the last case. Being real and symmetric, the function pi ((M|) has the 

Schmidt decomposition 

4 Asymptotic expansion of the RDM pi{x,y) ^^Xf'ui{x)ui{y), (27) 

Now we come to the main goal of this paper, which is 

to provide the Schmidt decomposition of the RDM in where {ui} and {AJ } are determined by 

the strong-correlation limit. To begin with our analysis, 

we take into consideration the RDM for the approximate / ,- \ r \j \('^) r \ fno\ 

^^ ' Pi[x,y)ui{y)dy := Xi 'ui[x), (28) 



wavefunction ijj'^ (|13|) . 

f .ap, \ lap, N, , /-.nN (u; |Mfc) = ^a- • One Can uotc that the introduction o 

Pap{x,y) = J ^B {x,X2,X3)i^B (y, X2 , 0^3)^X2^X3 (19) coordinates in the functions p2 (BSD and p^ m by 



which becomes exact as g —)• oo (xc ~> oo). An easy in- x i->- x — Xc, y H' y — Xc, (29) 

spection of Eq. ([T^ reveals that in this limit, it reduces 

to the form and 

X 1-^ X + Xc,y ^-^ y + Xc, (30) 

^"-P '^ Pi P2 P3, \ ) respectively, transforms them into g-independent forms 

^ff[f}i which are identical to each other, p2(a;,y),P3 (2;, y) 1— >■ p{x,y), 

Pi(x,y)=A/^.,(.r,x,Z2)^.,(zr,y,..)rf.id.., (21) p(x,y) « 0.2262e— (^^-^^^H™-. (31) 

The above function is also real and symmetric, thus its 
Schmidt decomposition is 



P2{x,y) = 

= A I 1pap{x + Xc,Zi,Z2)lpap{y + Xc,Zi,Z2)dzidz2,{22) 



p(i,y)-^Ap)i;z(x)«z(y), (32) 



;,(2) 
/=0 



where {vi} and {Aj } are determined by 



P3ix,y) = 

= A I i'ap{zi,Z2,X - Xc)lpap{zi,Z2,y - Xc)dzidZ2,{23) 



p{x,y)vi{y)dy = \f\i{x), (33) 

{vi\vk) = 5ik- By changing the variables back in ([5^ in 
where A = 2C^. We have obtained exact closed-form so- accordance with (EH), one gets the expansion of p2 in the 
lutions for the above integrals. However, since they are form 

quite lengthy, we report here only their numerical expres- p^(x, y) = V Ap'w;(x + Xc)vi{y + x^). (34) 

sions ;^o 

pi{x,y) « o.2407e"°-^^^'''''+°-^^^^^^"°-^^^''^'. (24) On the other hand, the change of variables back in (l32|) 

in accordance with pop . yields the expansion of p^ as 

PVsi^^y)- p^{x,y) = Y.^f\i{x-Xc)vi{y-Xc). (35) 

w 0.2262e"°-™^*^t(^'^^'=^^+(^='='^'=>^l+°-°'^^°(^^^''='(^='='''=\(25) '=0 

which does not hmit the generality of our further con- Obviously, the one-particle orbitals i;;(x±Xc) satisfy (wi(x± 

sideration. It it worthwhile to note here that the same ^^^(^ ^J^^t ^ ^ji:^^^^^^!^' substitution of Eqs. dUD, 

asymptotic behaviour for the RDM will be obtained for dM]), and dSS]) mto ^ gives, as g ^ 00 (x^ ^ 00), 
the case of fermions. g-i-oo _ 

Using Eq. (PH]) , we can obtain a closed- form asymptotic ^ ~ 

expression for the one-particle density: = > [Aj ui{x)ui{y) + A; vi{x + Xc)vi{y + Xc) + 

2 '=0 

n9^°° = 3p9^°°(x,x)«3(0.2407e-i-^389a; ^ 

^ \' ^ +A; V((x-Xc)ui(y-Xc)J. (36) 

-h0.2262e-i-^466(^-^=) -h 0.2262e-i-4'"^s(^+^'^) ). (26) 

In this limit, the family {ui(x),vi{x — Xc),vi(x + Xc)} 

An inspection of the computed 71^^°° shows that it ex- forms a complete and orthonormal set, since the integral 

hibits exactly the Gaussian peaks centred at the classical overlaps (^((x -I- Xc)|wfc(x — Xc)), {ui{x)\vk{x — Xc)) and 
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{ui{x)\vk{x + Xc)) vanish for any l,k. Wc can therefore 
recognize Eq. (p6| as the Schmidt decomposition of the 
asymptotic RDM p3^°° . Because the asymptotic natural 
orbitals vi{x + Xc) and vi{x — Xc) correspond to the same 

(2) 

occupancy A; , i.e., double degeneracies in the spectrum 
of the RDM p9^°° occur, the Schmidt decomposition (|36p 
fails to be unique [21] . Before going further we stress that 
the conservation of probability for the asymptotic occu- 



pancies gives X]z=o['^J + ^'^i ] ~ ^- ^'-'^ ^^^ ^^^^ °^ com- 
pleteness, we give below another form of the Schmidt ex- 
pansion of p^^°°, different from that of Eq. p6p . To be- 
gin with, we extend the results of [3T] to the case of more 
than one point of double degeneracy. From the orbitals 
Vi{z + Xc) and vi{z — Xc), which correspond to A; 
define the new orbitals to be 



(2) 



we 



and 



Vi{z) 



ri{z) 



Vl{z + Xc) + Vl{z~ Xc) 
Vl{z + Xc) -Vl{z- Xc) 

V2 



(37) 



that fulfill {'rii\Ti) = 0. In terms of them, Eq. ([55)) can be 
rewritten as 

p'^^ - 5][AWu,(x)uKy) + \^^^m{x)m{y) + 

+\f\i{x)n{y)], (39) 

Since in the limit as. g ^ oo {xc — >■ oo) we have (??i|%) = 
Sik-, {nW) = 5ik and the integral overlaps {ui\r]k): (w;|rfc), 
('7;|'''fe) vanish for any ^, fc, Eq. ((39)) yields nothing other 
than a Schmidt form different from Eq. ([M)) . 

The integral equations (|25)) and (|55)) can easily be 
solved through a discretization technique (see for example 
[IT]). The two lowest asymptotic occupancies are found 

numerically to be A^^^ « 0.3193, a[|^^ « 0.3249. Because 
the sum of all the remaining asymptotic occupancies. 



S/=i[-^/ + 2AJ '], is only about 0.03, it follows that in 
(|36p and in p9p only the terms with / = are important 
and, in consequence, p^~*°° approaches the form 

w X-^' UQ{x)uQ{y) + Xf'vQ^x -f Xc)vQ{y + x^) + 

+X\,'vf){x- Xc)vf){y-Xc), (40) 

in particular. 

Wc close our discussion with Fig. 01 which shows the 
numerically determined behaviour of A;, for I ~ Q — 2 
as a function of Ing. It is seen how the occupancies con- 
verge to their asymptotic values determined by the inte- 
gral equations (|28)) and (|33)) , which confirms their validity. 
In particular, we observe how Aq and Ai converge to an 
asymptotic doublet. 




Fig. 4. The dependence of the three lowest occupancies 
on Ing. The inset highlights their changes for large values 
of g. The asymptotic values are marked by horizontal lines. 
The values of asymptotic occupancies Ag ~ Ao^°° = 
^g-j-oo^ Ag = X2~^°° are reported in the text, whereas, the 
(38) values of three lowest occupancies determined as g — >■ 



with the use of 

0.6619, Ap° « 0.1755, Af 



16)) are found numerically to be Ag 
! 0.1004. 



in a ID harmonic trap. Using the harmonic approximation 
we explained the nature of the degeneracies appearing as 
g ^ oo in the spectrum of the RDM. An explicit represen- 
tation of the asymptotic natural orbitals and their occu- 
pancies has been derived in terms of some (7-indcpendent 
integral equations. Among other results, we found that in 
the g — >■ 00 limit the occupancies Ag ' = Aq^°° = Af ~*°° ss 

0.3249, X'-Q^ = Ag^°° « 0.3193 are the only three that have 
considerable values. In other words, it turned out that 
only the three natural orbitals contribute significantly to 
the asymptotic bosonic ground-state. We also determined 
numerically the three lowest occupancies as functions of g 
and showed how they tend to their asymptotic values. In 
particular, we obtained a closed-form asymptotic expres- 
sion for the one-particle density given as a linear combi- 
nation of Gaussian functions centred at the classical equi- 
librium positions of the particles. 

It would be interesting to fully investigate the effect 
of the number of particles on the correlation properties 
in ID systems of strongly interacting bosons and/or spin 
fermions with a Coulomb interaction. To the best of our 
knowledge, there is still a lack of studies along this line. 
We hope our results will stimulate others to undertake the 
investigation of this issue. This will be also a subject of 
our further research. 
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